COUPLING TECHNIQUES FOR NONLINEAR HYPERBOLIC EQUATIONS. I. 
SELF-SIMILAR DIFFUSION FOR THIN INTERFACES 
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Abstract. We investigate various analytical and numerical techniques for the coupling of non- 
linear hyperbolic systems and, in particular, we introduce here an augmented formulation which 
£\J ' allows for the modeling of the dynamics of interfaces between fluid flows. The main technical 

difficulty to be overcome lies in the possible resonance effect when wave speeds coincide and 
global hyperbolicity is lost. As a consequence, non-uniqueness of weak solutions is observed for 
the initial value problem which need to be supplemented with further admissibility conditions. 
This first paper is devoted to investigating these issues in the setting of self-similar vanishing 
viscosity approximations to the Riemann problem for general hyperbolic systems. Following 
earlier works by Joseph, LeFloch, and Tzavaras, we establish an existence theorem for the Rie- 
mann problem under fairly general structural assumptions on the nonlinear hyperbolic system 
and its regularization. Our main contribution consists of nonlinear wave interaction estimates 
, for solutions which apply to resonant wave patterns. 

-i— > 

1. Introduction 

This is the first part of a series devoted to analytical and numerical techniques relevant for the 
CNJ ■ coupling of nonlinear hyperbolic systems. We mainly discuss an augmented formulation which 

allows for the modeling of the dynamics of interfaces between fluid flows. The main technical 
difficulty overcome here for the Riemann problem (that is, a Cauchy problem with piecewise 
constant initial data) lies in the possible resonance effect when wave speeds coincide and global 
hyperbolicity is lost. As a consequence, non-uniqueness of weak solutions is observed for the 
initial value problem which need to be supplemented with further admissibility conditions. In 
the present paper, we restrict attention to self-similar vanishing viscosity approximations to the 
\ Riemann problem for general hyperbolic systems. 

Specifically, we are interested in the following class of nonlinear hyperbolic systems of (N + 1) 
partial differential equations 

Aq(u, v) dtu + Ai(u, v) d x u = 0, 
(L1) d t v = 0, 



O 



- 1—1 

X 



C3 



where the vector- valued field u = u(t, x) £ R N and the scalar function v = v(t, x) £ R (with x £ R 
and t > 0) are the main unknowns of the theory. We assume that the first set of equations in 
(|l.ip forms a strictly hyperbolic system but admits one wave speed that changes sign, so that the 
overall system is only weakly hyperbolic. Specifically, the mappings Aq, A\ are assumed to 
be smooth, matrix- valued maps such that Aq is invertible so that the first set of equations in (jl.ljl 
is formally equivalent to the following nonconservative system with variable coefficients: 

(1.2) d t u + Aoiu^v)- 1 A l {u 1 v)d x u = 0. 

It is assumed that the product matrix Aq(u, v)^ 1 A\{u, v) admits real and distinct eigenvalues 
denoted by Xi(u,v), 1 < i < N. Finally, it is assumed that one eigenvalue A m may take values 
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about the origin. For instance, there might exist a state (u*,v*) 6 x R such that the matrix 
A\(u*,v*) is non-invertible, with 

(1.3) X m (u*,v*) = Q. 

In view of the above assumption, the system f)l . 1 [) is called a weakly hyperbolic system. Our 
objective, precisely, is to study this resonant regime. 

Recall that Dafermos [IB [T3l E] advocated the use of self-similar regularizations in order to 
capture the whole wave fan structure of weak solutions to the Riemann problem. This consists in 
searching for self-similar solutions depending only on the variable £ := x/t and, then, introducing a 
self-similar regularization of the given hyperbolic system. Specifically, for the problem of coupling 
under consideration in this paper we propose, in the variable (x,t), to regularize (jl.ljl in the form 

Mu e ,v e )d t u e + Ai(u e ,» e )0 x u e = etd x (B (u £ ,v e )d x u e ), 

d t v* = e p td xx v e , 

where e > is a small parameter and Bq = B$(u, v) is a given matrix referred to as the viscosity 
matrix and p > is a real parameter. In the self-similar variable £, the equations satisfied by the 
viscous solutions (u e ,v e ) = (u e (£), v e (£)) read (with £ G R) 

( - £A)(w e ,v e ) +Ai(u c ,u e ))u| = e(B (u e ,v e )ul) e 

Our objective is to study the existence and regularity of solutions to (|1.4|) and to rigorously 
justify the passage to the limit e — > 0. Following earlier works by Joseph, LeFloch, and Tzavaras 
(see references below), we are going to establish a uniform (e-independent) bound on the total 
variation TV(u e ), and an existence theorem under fairly general structural assumptions on the 
hyperbolic system and its regularization. 

This general strategy was proposed in the case B = I and A = / by Tzavaras [34] (for 
conservative systems) and extended by LeFloch and Tzavaras [3D] (for non-conservative systems). 
The technique was further developed by Joseph and LeFloch in the series of papers [2B 1221 12H 
[2B [2B HZ] ■ For the purpose of the present paper, we will especially build on [27] where a general 
technique to derive interaction estimates was introduced and general matrices Bq were dealt with. 
For other results on self-similar limits including viscosity-capillarity terms and large data, we refer 
to pioneering works by Slemrod [16], Slemrod and Tzavaras [33], Fan and Slemrod [15] . 

The coupling that we are studying in the present work may be non-conservative in nature 
(cf. Section [2] for details). In contrast, for coupling techniques based on systems in conservation 
form, a large literature is available; see for instance [3 IB IB IB fTH I3T] , 

An outline of this paper follows. In Section [B we present a general approach involving the cou- 
pling of nonlinear hyperbolic systems. In Section [3J we discuss the case of scalar- valued unknowns 
u, and establish a general existence theorem for the viscous self-similar Riemann problem (11.41) 
when N = 1. As e — > 0, we prove that this smooth solution converges to an entropy solution, 
at least in each half-space x < and x > 0. This global existence result requires no smallness 
assumption on the data nor on the coupling of the two models. 

The core part of this paper is contained in Sections [4] and [5] which cover general systems of N 
equations. Imposing a natural smallness condition on the Riemann data and the coupling of the 
two models, we establish the existence of smooth solutions to the viscous problem, even in the 
presence of a resonance effect. In Section @] we derive the main estimates on the total variation 
while, in Section [B we justify the limit e — > 0. We refer to the forthcoming works [B IB 13 f° r 
further investigation of these solutions and related issues. 

2. A FORMULATION OF THE COUPLING OF NONLINEAR HYPERBOLIC SYSTEMS 

Before we can state our new formulation based on an augmented system, we start by briefly 
explaining the formulation based on a fixed interface. The weakly hyperbolic problem mentioned 
in the introduction arises, in particular, via the following coupling technique. Consider two strictly 
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hyperbolic systems posed in half-spaces: 

d t w + d x f-(w) = 0, x<0, t>0, 
d t w + d x f+{w) = 0, x > 0, t>0, 

where the flux f± are given smooth maps defined on open subsets f2± C R and the unknown of 
the problem is w — w(t, x) G fi_Uf2+. In addition to initial data, a certain coupling condition must 
be prescribed at the (fixed) interface {x — 0}. This problem can be regarded as a boundary and 
initial value problem within each half-problem, and the fundamental question is how to formulate 
a physically relevant boundary condition so that the global problem is well-posed. One natural 
requirement, following Godlcwski and Raviart [ElllB], is imposing the continuity condition 

(2.2) 9-(w(0-,t))=9+(w(0+,t)), f>0, 

where , 6 + are two invertible maps in I 1 *' , with inverses 

These functions precisely provide the necessary freedom to handle various types of couplings. For 
example, by choosing 6- = 6 + = Id one imposes the continuity of the variable w at the interface, 
while by choosing 9± = f± one imposes the continuity of the flux at the interface (so that the 
general problem is conservative). 

Recalling Dubois and LeFloch's theory [15] of the initial and boundary value problem for 
nonlinear hyperbolic systems, it is clear that the condition (|2.2[) is realistic only when the boundary 
is not characteristic, that is when all eigenvalues are bounded away from 0. In the latter case, 
instead, following the weak formulation of the boundary conditions proposed in |15) generalized, 
for the coupling problem, by Godlewski and Raviart [TH HH] and Ambroso et al. [2], we impose 
that the interface condition is satisfied in a weak sense, only and, specifically, takes the form 

w{0-,t) G $_(0_o0+ 1 (w(O+,t))). 

where $ + (6 + ) (resp. $_(&_)) is the Dubois-LeFloch's set of admissible traces of the associated 
Riemann solutions 

$+(M := {i?+(0+;fo+,a), a G 0+}, 
$_(&_):= {i?_(0-;a,&_), a G Q-}. 
Here, R = R+(x/t; &+, a) denotes the solution of the Riemann problem 

d t R + d x f+{R) = 0, xel, i>0, R(x,0) = 

and similarly R = R_(x/t; a, &_) is the solution of 

d t R + d x f_(R) = 0, i>0, R(x,0) = 

Yet, when /_ =/= f+, the question of the existence and uniqueness of weak solutions satisfying 
(|2.3p is a challenging issue. In the present work, we propose to reformulate the above problem by 
"removing" the interface and defining a new problem posed on the whole space R. 

We proceed as follows. First of all, we define the new variables 





(2.4) U-:=9-(w), u + :=0 + (w), 



I U-, x < 0, 

u := < 

I u+, x > 0, 

and we rewrite the half-space problems in the (conservative) form 

(2.5) d a ±{u) + d x f±{j±{uj) =0, ±x > 0, t > 0, 
or equivalently in the (nonconservative) form 

(2.6) {D ul ±{u))d t u + {D 1 f±){ 1 ±{u)){D ul± {u))d x u = Q 1 ±x > 0, t > 0. 
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The coupling condition becomes 

u(o+,t) e*+(u(0-t)), 
(2 ' 7) «(o-,t) e *_(«(o+,t)), 

where ^+(b) (and similarly is the following set of admissible trace at £ = 0- 

*+(&) = {i?+(o+,&,u + ), u+ g n}, 

and it = 6, u+) is the self-similar solution of the following Cauchy problem 

9 t7+ (i t ) + 9 x / + ( 7+ (w)) = 0, i6M,(>0, u(a;,0) = 

In absence of a resonance phenomenon, this reformulation allows us to simply impose the continuity 
of u at interface 

(2.8) u(0-,t)=u(0+,t). 

Second, we propose to replace the problem (|2.5p - (|2.7|l by the new problem (already mentioned 
in the introduction) 

Aq (u, v)dtu + A\ (u, v)d x u = 0, 
(2 - 9) d t v = 0, 

which is a nonlinear hyperbolic system in nonconservative form |281 129| and where v : [0, +oo) x 
R — > [—1,1] will be referred to as the color function. We arrange that regions where v = — 1 
correspond to the left-hand half-problem while regions where v = 1 correspond to the right-hand 
half-problem, by requiring the following consistency property on Aq, A\: 

A (u,±l) = D u ~f±(u), 

(2 ' 10) A 1 (u,±l) = D 7 / ± (7 ± («))£> u 7±(«). 

and assuming the existence of a function C — C{u, v) so that 

A (u,v) = D u C(u,v), 

and C(u, ±1) = j±(u). For j = 0, 1, by definition, the matrices Aj(u, v) should smoothly connect 
Aj(u, — 1) to Aj(u, 1) as v describes the interval [—1,1]. Moreover Aq must be invertible and 
Aq 1 A\ have real and distinct eigenvalues, extending here the strict hyperbolicity of the original 
hyperbolic half-problems. 

The system (|2.9p is then supplemented with the initial data 



,2n\ u o (x,0) =u (x) =: 9±(w (x)), ±x > 

Vq(x, 0) = Vq(x) := ±1, ±x > 0, 

for some given data uq. 

We are especially interested in the case that the interface is characteristic for some state value 
(u*, v*), that is, when the matrix A\(u*, v*) admits the eigenvalue and (|2.8|) need not be satisfied 
as an equality, in general, so that the weak formulation above is necessary. 

3. Existence theory for scalar conservation laws 

3.1. Riemann problem with diffusion. In the present section, we restrict attention to scalar 
equations and provide a rather complete study of the problem described in the introduction. Note 
that the problem under consideration is nonconservative in nature, and reduces to a conservative 
system when the component v takes the values ±1. As explained earlier, we search for a function 
U = w(£) obtained as the limit of smooth approximations u e ,v e to 

{-ZM(u € ,v*)+A- L {u*,v e ))u\ = e (Bo(u e ,« e )«!) e , 
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supplemented with Riemann initial data 

u £ (— oo) = ml, u c (+oo) = ur, 

(3 2) 

y ' 7 w £ (-oo) = -1, v £ (+oo) = 1. 

In (|3.1[) , the maps Ao and Ai are now smooth scalar-valued functions, which satisfy the follow- 
ing consistency condition with the underlying hyperbolic coupling problem determined by the 
functions 7 ± and f± : there exist constants c\ , C2 , C3 , such that 

(3.3) < ci < i4o(«, «), < c 2 < S (m, w) < c 3 , 

and 

A (u, -1) =7-H, 4>(u,l) =7U U )> 

i4 1 (tt,-l) = (/_o 7 _)'(«) 1 A 1 (t»,l) = (/+o 7+ ) / (ti). 

We set £1 := [mm(uL, u.r), max(ii£, Mr)], and introduce the Lipschitz constants uo, u>i of A$,A\, 
respectively, i.e. 

v) — Aj(u, v)\ < uij (\u — u\ + \v — v\) 

for all (u, v), (u, v) £ X [—1, 1] and j = 0, 1. 

The first equation in (|3.1|) can be equivalently rewritten as 

(3.5) (-£ + A(<<))GV,<>A)(^<H4 = e {B Q (u\ , 

where A and G are defined by 

A {u,v) B (u,v) 

Furthermore, our assumptions imply that (for some A > 0) 

|A(u,v)| < A, (it, v) GOx [-1,1], 

which expresses the property of finite speed of propagation for the underlying hyperbolic problem. 

Given any M > A, we will study first the problem on the bounded interval [-M, M] with the 
boundary conditions in (|3 . 2[) imposed at the end points ±M. Later, we will let M tend to infinity. 



Proposition 3.1 (Existence for Riemann problem with diffusion). For each e > the problem 
(|3.ip admits a smooth solution (u £ , v e ) £ C°([—M, M], fl x [—1, 1]) (space of continuous functions) 
given by the implicit formula: 



(3.6) 



e -fc'( u «iO/« B (u e ,v e )- 1 dC 

-M 

M ' 

e -fc'(«« 5 0/ e B {u € ,v e )- l dC, 

—M 



m £ (£) = -1 + 2 



? e-^/ 2£P dC 



— OO 



+ OO 

e^ 2 / 2£? d( 



with 



h e (u e ; £) := / (C - A (m £ , « £ ) ) G( U £ , « £ ) 



Moreover, these solutions u e and v e are monotone, bounded, and continuous, and have uniformly 
bounded total variation: 

TV(u e ) <\u R -u L \, TV> £ )<2. 
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Proof. Solving the second equation in (|3.1[) is immediate. On the other hand, we can rewrite the 
problem (|3.ip as 

B Q (u c ,v e )u\ = ip, 
(-£ + A(u> 6 )) G(u e ,v e ) tp = e<pt. 

Given u G C°([— M, M], 0) we consider the solution u e (u;£) of the "linearized" problem 

Bq(u, v e )ul = ip, 
(-Z + \(u,v e )) G(u,v e ) (p = etp c , 

together with the boundary conditions 

u e (-M)=u L , u € (M) = u R . 

The solution is explicitely given by 

f ^(«;C)Bo(«,« e ) _1 dC 

J-M 



u £ (u; f) = u L + (u R - u L ) 



M ' 

, p{u-OB (u,v f -)- l dQ 

{6.1} J-M 

<p(u;£) = exp ( - /i e (u;£)/e), 
h e (u;0= f\<;~X(u,v e ))G(u,v t )d(, 

J a 

in which a G [— M, M] is arbitrary. The above formulas determine a map T £ that takes u G 
C°([-M, M], fi) to the function u e \u\ •) G C°([-M, M],Q). We need to hnd a fixed point of T e . 

The uniform bounds on A (it, u) and G(u, v) (for any u G C°([—M, M], and v6 [-1,1]) allow 
us to choose a e G [— M , M] so that 

h e (u;^)>0, £e[-M,M]; h t {u- ) a t )=0. 

Consequently, for all £ G [— M, M] and for some constant c 4 we have 

0<h e (u;O <c 4 , 

so that 

C3 1 exp (-c 4 /e) < ^(G^SqCM 6 ) -1 ^ < c^ 1 . 
We also obtain the uniform bound 



< \u R - u L < 



M 

-M 

-1 



< \Ur — Ul\ ~i 

2Mc 3 exp(— c 4 /e) 



< \ur " u L 



"3 

c 3 exp(c 4 /e) 
2Mc 2 



The bound above being independent of u, we deduce that the family T e is equicontinuous and 
its image is relatively compact in (7°([— M, M],U\. Since this image is a convex closed subset 
of the Banach space C°([— M, M]), Schauder's fixed point theorem applies and ensures that T e 
admits a fixed point. Hence, there exists u e G C°([— M, M],f2J such that T e (u e ) = u e , and 
the representation formula (|3.6p holds. The uniform total variation bounds follow directly from 
(EU). □ 
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3.2. Passage to the limit. Using the notation in Proposition ^. 11 we continue with the following 
two lemmas. 

Lemma 3.2 (Existence of a pointwise limit). After extracting a subsequence if necessary, the 
sequence u e converges pointwise to a limiting function u lying in the space BV([— M, M]) (of all 
functions with bounded variation): 

u e (£)^u(0, £e[-M,M], 

which satisfies, in the sense of distributions, 



(3.8) 



•£^7-(«) + ^/-(7-(«))=0, £<o, 
•£^7+(«) + ^/ + (7+(tO) = 0, e>0. 



Lemma 3.3 (Entropy inequalities). The limit u given by Lemma \3.2\ also satisfies, in the sense 
of distributions, 



(3.9) 



-e^»?(7-(«)) + ^g-(7-(«))<0, £<0, 

-e 4»?(7+(«)) + 4«+(7+(«)) < o, e>o, 



for all convex entropy functions r\ and associated entropy flux q' ± = rf f. 

Proof of Lemma lKM In view of Lemma \3. 11 Helly's compactness theorem applies and, as e — > 0, 
ensures the existence of a pointwise limit (u,v) with bounded variation. Fix 9in(0,M) and let 
4> G Cq°((9, M)) be a compactly supported test-function. In the integral form, (|3.1|) becomes 



A I 



A I 



£ A ( U »| </>d£+ / Ai(u c ,v c )u\4>dt = e (S (u»f)<M£, 



A I 



that is 



where 



A I 



At 



£ A (u e ,l)u\ </>d£ + / Ai(u e ,l)«| ^df + fi e = e / (B (it e ,« e )ttf) f 0d£ 



A/ 



A/ 



$ 0W,1)- A,( U > e )) U | ^d£- / (Ai(u e ,l) - Ai(u e ,?; e ))u| 0d£. 



A/ 



Using (|3.3p and (|3.4p . we can write 



A/ 



A/ 



/o ^ 

The term Q e vanishes with e, since 

i-M 



/+( 7+ K)) 4>dt + W = e (B (u e ,« e )u|) 0d£ 



A/ 



|ft e |<M / w |l -« e | |0|^+ / wi|l - v e \ |</>|^ 
<(Mw +Wi) |l-?/(#)| II^HooTU^), 

where the total variation term TF(u e ) remains bounded and |1 — w e (6')| tends to 0. On the other 
hand, we have 



A I 



{B (u e ,v e )ul)^dt; 



A I 



(B (u e ,v e )u\) c de 



<e ll^lloo c 3 TU(u £ ), 
which also converges to zero. As e tends to 0, we conclude that 

cA I 



J r-e^7+(«) + ^/+(7+(«)))0de = o, 
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which is the first condition in (I3.8[) . The same arguments apply on the interval [-M, 0], by using 
test-functions supported in the interval (— M, 9), with 6 < 0. □ 

Proof of Lemma VJJA Fix 8 > and let <j> £ Cg°((0, M)) be a non-negative test-function. Multi- 
plying (|3.ip by ?7'(C(u £ , 1))0, we get 



(3.10) 

Observing that 



£A (ti%<)K|j/(C(</,l)) + / -4i(u e ,« e )u|?7 (C(u £ ,l)) <M€ 



A I 



(B (u*,v*)ul) rf (C«l)) tfd£. 



£ 0W, « £ ) - A (i/ e , 1)) ^'(C(« e ,l)) 0d£ 



<Mwb|l-i; e (e)| TV{u<) U{C{-,1) 



I oo I V 7 II oo i 



and similarly for the coefficient A\, we see that the left-hand side of (|3.10[) is equivalent (modulo 
terms that tend to zero with e) to 

£d u (C(u e , 1)) ulr]'{C(u £ , 1)) 4>di + / A!( u e , 1) uf»/(C(u e , 1)) 0d£ 

Jo 

w d , . d 



On the other hand, the right-hand side of A3. 10)) can be rewritten in the form 



M 



(B (u\v')ul)^'(C(u',l)) <pd£, 

M 

B a (u e , v e ) (ul) 2 A Q (u\ l)v"(C{u e , 1)) 0d£ 

M 

B Q {u\v e )u\q'{C{u e ,l)) <^d£, 

o 

in which the first term is non-positive and the second one tends to 0. Thus, letting e —> we 
obtain 

rM d r M d 

^- V ( 1+ (u)) <t><% + J q ^g+(7+(«)) <t>*Z < o, 

which yields the second identity in the statement of the lemma. The derivation of the first identity 
in the half-space £ < is completely similar. □ 

3.3. Riemann problem for the hyperbolic coupling problem. In view of the boundary 
condition (13.21) . it is natural to extend u by 

V, £<-M, 
u R , £>M. 



(3.11) u(0 



The conclusions in Lemmas l3.2l and r3.3l then clearly hold on the intervals (-co, —M) and (M, +00). 
In addition, an interface condition for the solution u at the end points £ = — M and £ = M is now 
derived, which is necessary to ensure that (I3.8[) and (|3.9|) extend (in the sense of distributions) to 
(0, +00) and (—00,0). 

Lemma 3.4. As e tends to 0, the solution u e converges uniformly toward ur (repectively ul) on 
the interval (A, M] (resp. (-M,—A[). 
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Proof. Let f G (^±^,M) be given. According to (|3.6I) we have 



K(0 - ur\ = \ul - u R \ 



< \ul - u R \ 



M 



M 



-M 
M 



i- he ^' e B {u e ,v e )- x d£ 

e -h<(u*)/e Bo ( u e jV eyl d £ 

e - h '( u °y e B Q {u e ,v e )- l d£, 



M 



M 



We use here the constant a G [—A, A] as the lower integration bound for both h e and the function 

HO ■= [ (C - A(t»(0, «(C))) G(u(C), w(C)) dC, 



hence h > 0. Moreover, h e (u e , ■) converges uniformly to ft, with 

|/i e (u e ,0-M0| 



(A(« e (0, « e (0) - A(«(C), «(C))) GK(C), v'(0) d( 

< + Pilloo^) ||G|U J" (|« e (C) - "(01 + |» e (C) - KOI) dC 

< Q«l + Pllloo^) ||G||oo(|K - U \\ L l + \\v< - v\\ L1 ). 

The uniform convergence of /i e (u c , •) toward a positive continuous function h such that h(a) = 
insures that there exists eo > together with A > B > and r\ > such that, for all e < cq, 

h e (u e ,£)>A, A+ n M <£< M, 



Thus, we deduce that 



h e (u c ,0<B, £ G [— M, M], |£-a|<>j. 



M-A c -A/e c -l 

|u c (0 - < - wl| 2 



,M-Ac 3 _, A _ 



2f) c 2 



(A-B)/e 



and u e converges uniformly toward u r on the open interval (^f-.M). The same argument leads 
to the uniform convergence of u e toward mj, on (— M, A ~ M ) . 

We summarize the results established in Lemmas 13.2113.31 and 13.41 as follows. 



□ 



Theorem 3.5 (The Riemann problem for the coupling of two scalar equations). Up to extracting 
a subsequence, the solution u e to ()3.1j) - (|3.2j) converges pointwise to a function u G BV(M), 

which satisfies the conservation laws and entropy inequalities 

-e^ 7 -(«) + ^/-(7-(«)) = o, e<o, 

-S^7+(t») + ^/+(7+(«)) = ) e>o, 



in 
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and 

-£^*7(7-(«)) + ^g-(7-(«))<o, £<o, 

-£ ^( 7 +(m)) + ^g+( 7 +(«)) < 0, £>0, 
/or aZZ convex entropy pairs, together with the boundary conditions 

u(— oo ) = ul, it (+oo) = ur. 
Equivalcntly, in terms of the function iu in (|2.4[) , we have established 

-f^« + ^/+(«0 = o, -e^H + ^y+H <o, £>o, 

with 

io(— oo) = wl, w(+oo) — wr. 
We have thus established that the interface problem admits a solution which has bounded variation. 

4. Existence theory for systems 

4.1. Terminology and notation. We will now generalize the technique developed in Tzavaras [31] 
and Joseph and LeFloch [22, 25], and cover the class of nonconservative and resonant systems p. 41) 
under consideration. We follow closely the notation and presentation in [25 . 

Specifically we consider the diffusive Riemann problem (|1.4I) with Riemann data ul,ur, and 
establish that, provided ul, ur £ 23(o"i) with a sufficiently small 8i < Sq and under some structural 
hypotheses on the matrix fields Aq,A\, this problem admits a smooth, self-similar solution, u e — 
u e (x/t) £ 23(00) and v e = v e (x/t) £ [—1,1], which has uniformly bounded total variation 

(4.1) TV(u e )+TV(v e )<C, 

for some uniform C > 0. Solutions to (|1.4[) will be decomposed in terms of "wave strengths" of 
the associated Riemann problem (|1.1[) . The uniform estimate (|4.1[) is the key to the convergence 
analysis as e — > 0, and the proof of the existence of the Riemann solution to the underlying 
hyperbolic problem, discussed in the following section. 

We are interested in solutions u taking values in a small neighborhood of a given state (normal- 
ized to be the origin without loss of generality), that is, in the ball 51 := 23(o"o) with (small) radius 
Sq. For each u £ £(oo) and v £ [—1,1], let Xi(u,v) < ... < Xpr(u,v) be the real and distinct 
eigenvalues of the N x N matrix 

A(u,v) := Ai(u,v) Ao(«,u) _1 , 

and let l\{u, v), . . . , In( u i v) and r\(u, v), . . . , r^{u, v) be basis of left- and right-eigenvectors, re- 
spectively, normalized so that U(u,v) ■ Vj[u,v) = if i 5^ j and li(u, v) ■ Ti(u,v) = 1. 

By reducing Sq if necessary, we may assume that the wave speeds A,-(tt, v) are sufficiently close 
to the constants Ai(0,0) and, in particular, are uniformly separated for all u £ "B(So) in the sense 
that, for some constants 

—M < A ± < Ai < Aa < . . . < Ajv < A N < M, 

(4.2) Aj := A,(0, 0) - O(S ), A,- := A,(0, 0) + O(S ) 
and 

(4.3) < Xi(u,v) < A,, tt€2(*o) s »€[-l,l]. 

Let m be the index associated with the resonant wave, i.e. such that X m may change sign. In 
addition, for So sufficiently small the vectors ri(u, v) are sufficiently close to n(Q, 0) and we assume 
that 

k(ui,v) ■ n(u 2 ,v) > I ~ S , 1*1,112 £ 2(o~o), v £ [—1,1], 

(4.4) 

|Z,-(«i,u) -rj(u2,v)\ < Sq, ui,u 2 £ T>(8o), i^ j, v £ [-1,1]. 
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In (|1.4p . the matrix B = B (u, v) is assumed to be non-degenerate and depend smoothly upon 
u and v. We treat the case that the diffusion matrix B(u,v) := Bq(u, v) Aq(u, v)^ 1 is sufficient 
close to the identity matrix, that is, for some given matrix norm and 7? > sufficiently small, 

(4.5) sup \B(u,v)-I\<Tj. 

ues(<5 ) 

«e[-i,i] 

To handle arbitrary diffusion matrices B(u,v), we follow Joseph and LeFloch |27] and introduce 
the generalized eigenvalue problem: 

(-£,ld + A(u,v))r t (u,v,£,) = m(u,v,^)B(u,v)ri(u,v,^), 

(4.6) ^ 

k(u, v, £ ) • (-£ Id + A(u, v)) = fafa, v, £) k(u, v, £) • B(u, v). 

with unknowns the vectors r»(u, v, £), h(u, v, £) and the scalars fii(u, v, £). We impose the following 
normalization to generalized left- and right-eigenvectors: 

n(u,v,£) ■r i (u,v,^) = 1, 
li(u,v,g) ■ B(u,v)rj(u,v,£) = if * ^ j, 
k(u, v, £) • B(u, v) ?i(u, v, f) = 1. 
Multiplying the first equation in (|4.6p on the left by r~i(u, v, £) and rearranging terms, we get 

(4.7) /^(u, v, = (-£ + Aj(u, u, 0) V) 0) 
where 

^ 4 ^ Ai(u,u,£) := ri(u,v,£) ■ A(u,v)?i(u,v,£), 

l/di(u,v,£) := fi(u,w,f) • B(u,v)?i(u,v,^). 
Clearly, in the special case where B(u, v) = /, we find 

Hi(u,v,£) = -f + Xi{u,v), ?i{u, i>,£) = n(u,v), k{u,v^) = k(u,v). 

So, by continuity, when B gets closer to the identity matrix, the coefficients dj(it, and 
Ai(u, v, £) get closer to 1 and Xi(u, v), respectively. In consequence, under the assumption \B(u, v)~ 
I\ < r\ with rj sufficiently small and by increasing the gaps Aj — A^ if necessary, we can always 
assume that 

h-Oiv) <%{u,v,C)<^i + 0{r,), 
l-0{rj)<di{u,v,() < 1 + 0(7?) 
for u £ 23(<5 ), v £ [-1, 1], £ £ [-M, M]. The following property was pointed out in [27] . 

Lemma 4.1. The ^-derivatives of the generalized eigenvectors and eigenvalues satisfy 

1^(^,01 = 0(7?), 

( ' ' dtiH(u ) v ) t) = -l + 0(r } ). 

Finally, we introduce the coefficient 

(4.10) v :=sup h -d v (Bfj) 

allows one to measure how closed the left-hand and right-hand hyperbolic models are and, from 
now on, this coefficient is assumed to be sufficiently small. 

Example 4.2 (p-system). The following example illustrates the meaning of v. Consider the coupling 
between two systems of two conservation laws, specifically p-systems, with two different pressure 
laws p± = p± (t) 

d t r - d x V = 0, 
d t V + d x p±(r) = 0, r>0, V£R, 

in which the associated flux F± have Jacobian matrices 
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Suppose that the coupling is based on B — Id and the average matrix 

Then, the eigenvalues and eigenvectors of the system are (j — 1,2 and ± corresponding to the two 
wave families) 



A,.r.r. = ^ -I±V + ( T )_ ^yV_(r), 



f i (r, U )=(l, T ^_i±2£ p ;( T ) _ 

/>,«) = (±^_1±V + (t) - l^p'_(r), l) , 
and a tedious calculation yields 

\p'(r)-p'_(r)\ 



h ■ d v (B?j 



< 



< 



2 s J- 1 -^p / + (r)- 1 -^p'_(r) 
\p' + (r) - p'_(r)\ 



2 min U/-p+(r), V / Z p r (r)') 
so that 

v < C sup \p' + —p_\. 

4.2. Equations satisfied by the characteristic coefficients. We supplement (cf. fll.4[) ) 

^(« e ,« e ) + Ai(« e ,« e ))«| - e(£ (u»f) s , 



(4.11 



with the following boundary conditions, inherited from Riemann initial data 

u e (-M) = u L , u e {M)=u R , 
*~ ' ' v e {-M) = -1, v c {M) = 1. 

We describe now an asymptotic expansion for the solution of the diffusive Riemann problem (|4.1ip - 
(|4.12[) . To handle an arbitrary diffusion matrix, the decomposition must be based on the modified 
eigenvectors (|4.6[) . We first solve explicitly the equation concerning v e : 



_ s J 

(4-13) ^(fl:="f(fl = 2 " ^ 2 ■ 

J_ M e ^dx 

Remark 4.3. In the limiting case p — +oo, we would formally get a Dirac mass solution if) = = 
2<5^ = o and the transition from vl = — 1 to vr = 1 at t/ie interface would then be discontinous. 

The main equation (|4. 1 1[) now reads 

(4.14) (-£T + A(n £ , A Q (u\ « £ )u| = e(i?( M e , i» e )A ( U e , «>f) { . 

To deal with this equation, we introduce a decomposition of the vector Ao(u e , v c )u|(£) on the 
basis of eigenvectors rj(u*(£), i> £ (£), £), that is, 

(4.15) ^ 
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where the functions a* are referred to as the characteristic coefficients. Removing the explicit 
dependence in e, the right-hand side of (|4.14p takes the form 

(B(u, v) Ao(u, v)u i ) ( = ^2 a' k B(u, v) ?j{u, v, ■) 
j 

+ ^ a j a k Du(B rj)(u, v, ■) A (u, v)' 1 f fe (u, v, ■) 

+ ^2 a J d v{Brj)(u, v, ■) V£ + 2J aj B(u, v) d(fj(u, v, •)■ 
j j 
Therefore, given any solution to (|4.11[) . we obtain 

^2 ( ea '' j B(u,v)rj(u,v,-) ~aj + A(u,v))rj(u,v, •)) 
j 

= -e ^ ajak D u (B?j){u,v,-)A (u,v)~ 1 r k (u,v>-) - Oj d v (B r )(u, v,-)ip 

3,k j 

-e^dj B{u,v)d i r {u,v, ■). 

3 

Now, multiplying the above equations by each vector li(u, v, •) for i — 1, . . . , N and relying on 
the equation (|4.6[) . we arrive at a coupled system of N differential equations for the characteristic 
coefficients cii'. 

(4.16a) a- — — — a, = rjL^u, v, •) + Qi{u, v, •) + S t (u, v, ■), 

e 

where the linear, quadratic, and source terms are defined by 

Li(u,v, •) := y^^jj(u,v, -)a,j, Qi(u,v, •) := ^ K ijk(u,v, ■) a a k , 

(4.16b) ' 

Si(u,v,-) := 2_^(Tij(u,v, -)ajip, 
j 

respectively, with 

Wij(u,v, •) := -r)~ l (h ■ Bd^rj)(u,v, •), 

(4.16c) Kijk(u,v,-) := -(k ■ D u (Brj) A 1 f fc )(w,u, ■), 

<Tij(u,v, •) := (h ■ d v (B?j))(u,v, •)■ 

The main equation (j4. 1 1[) is therefore completely equivalent to f|4. 1 5[) - (|4. 16[) . In view of (14.91) 
and (|4~TU|) , (t4~l%|) takes the form 

(4.17) a' i --^ i (u,v,-)a i = 0(r 1 ) + O(l) £ | 0i | |o fc | + 0(i/) £>,-| |# 

3 j-k j 

Consider the principal part of (I4.16[) , that is, given some function u = u(y) (which at this stage 
need not to be a solution to (I4.16[l ) let us consider the following decoupled system of N linear 
equations: 

(4.18) ^_/^ii^ = , i = l,..., N . 

The general solution is (up to a multiplicative constant) 

e -ffi/« rv 

(4.19) if* := , g t (y) ■= ~ I (M(u,v,-)(x)dx, 

/ e~^ e dy Jpi 

J-M 

where the constants pi will be chosen so that the functions gi are non-negative (cf. Section 
Clearly, the functions ip* are strictly positive and their integral over R equals 1. 
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We will search for the general solutions a t of (|4.16a[) in the form 
(4.20) a z = Tttf + 0i 

where Tj refers to a wave strengh and with 9i small w.r.t n. When solving the equation (I4.16a[) for 
a given right-hand side and considering again first-order terms, we are naturally led to consider 
the following linear wave coefficients Jj->i, the quadratic wave coefficients Fjk^i and the resonant 
quadratic coefficients Jj_^ defined by 

(4-21) Jj^M := <pt(v) / ^ri ^ 

Jci Pi \ x ) 

(4.22) f^m-vUv) f 

(4-23) Jf^(y) := (£) [ dx 

for some constants c, € [A i; Aj] independent of e. By studying these coefficients, we will gain useful 
information on the possible growth of the total variation of solutions to (|4. 16|) : roughly speaking, 
Jj-yi bounds the influence of the j-th family on the i-th family, Fj_k-*i bounds the contribution 
on the z-th family due to waves of the j-th and A:-th characteristic families and bounds the 
influence of the j-th family on the i-th family "through" the coupling wave if). 



4.3. Linearized wave measures. The results of this section were established earlier in Joseph 
and LeFloch [22] and are presented for the convenience of the reader. We fix some speed range 
[A m in, A max ] and analyze the formula (|4.19j) within this speed range. We introduce first the space 
of functions that are "almost linear at infinity". Observe the coefficients //, will belong to such 
spaces in their respective speed ranges [A,-, AJ. 

Definition 4.4. A function h : [-M, M] i— > R is said to belong to the class L of almost linear 
functions if there exists two functions d, X G L°°([— M, M],M), and two positive reals <i m m, <inax 
such that 

h(x) = d(x)(X(x) - x), xe[-M,M] 

and 

< e?min < d(x) < 

"max i 

-M < A min < X(x) < A max < M. 
Lemma 4.5. Let h : [— M, M] H y M be a function of class L and, given y G [— M, M], set 

g(x) = - [ h(x')dx', xe[-M,M]. 



J y 

Then, g is Lipschitz continuous and achieves its global minimum at some (non-unique) point 
9 S [A min , A max ] such that X(p) = p. 

Proof. Using definition 14.41 of the class £ we get 

h(x) > 0, X < A min , 

h(x) < 0, X > Amax, 

g being continuous, decreasing on [— A/, A m i n ] and increasing on [X max ,M] (since g' = —h), we 
deduce g achieves its global minimum at some (non- unique) point p £ [A m ; n ,A max ]. Moreover 
h(p) = —g' = 0, that means X(p) = p. □ 
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We now introduce a new notation suggested by the formula (|4.19l) and derive useful algebraic 
properties. For x,y G [-M,M] and ft £ £*([— M, M], K), we set 

4>(y, x; ft) :— exp ^— J h(x') dx'^j , 

I{V,h):= f <p(y,x';h) dx', ip(x; ft) := ^ . 

That is, x h4 <^(a;; ft.) is a solution (with unit mass) to the differential equation 

, ft 
ip' - ~ip = 0. 
e 

Lemma 4.6 (Algebraic properties of the mapping <fi). For general functions h,h in i 1 ([— M, M], R) 
and x,y, z G [— M, M] £fte following algebraic properties hold: 

4>{x,y;h) - 

— =7- = cj>(x,y;h- ft), 

4>{x,y;h) 

(f>(x,y;h) <p(y,z;h) 



>(z,y;h) 4>{y,x\h) 
(x,y;h) = <f>(y,x;-h) 



x, z; ft), 



<j)(y,x;h) 

In view of these properties, the definition of y>(a;; ft) is checked to be independent of the variable 
y. When the function ft belongs to the class L, by taking y — p given by Lemma 14.51 we obtain a 
negative argument in the exponential defining (f> and that argument vanishes at points where g is 
minimized. 

Now, fix /xi,/X2 in the class L, and let tp\ := ip(-\ /ii) and if2 '■= </?(•; ^2) be solutions to the 
differential equation associated with p\ and fi2, respectively. Denote by p\ and P2 the minimization 
points of the associated functions as defined in Lemma 14.51 Moreover, fix some yet unspecified 
scalar c G [A m ; n , A max ] which we take to be independent of e. We then want to control the linear 
coefficient 

(4.24a) J^viv) = <Pl(v) f dx 

characterising the first order (linear) influence of ip2 on ip±. Using Lemma 14.61 we find those two 
following equivalent expressions, both useful in the sequel according to the sign of \i% — 

J<p 2 -»pi{y) = 77^-4 ^(p2,pi;Ma)yi(y) / <i>(pi,x;n2 - m) dx, 



J<P2^<pAy) = ¥2{y) J 4>{x,y,m- H2) dx. 

In order to estimate those coefficients, we also need for more information on the asymptotical 
behavior of appearing quantities as e tends to 0. The two following lemmas will give it. 

Lemma 4.7 (Asymptotic behavior of (j)). Let [x,y] be an interval of '[—M, M], (withx < y), andlet 
ft be a continuous function on [x,y]. Ifh is strictly positive bounded, say we have h(x') > ft m in > 
on [x,y], then the following integral is at most linear in e 

v e 



(x , y, —ft) dx' < 

Lemma 4.8 (Asymptotic behavior of I(p; ft)). For a function h in the class L, and p defined in 
Lemma \4^j the integrals I{p;h) satisfy 

ce < J(p; ft) < 2M. 
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Proof of Lemma \4- 7j Let us suppose h(x') > h min > 0. Then, for all x < x' < y, we have 

h(t) dt < -hmin (y - x'), 



thus 

j>(x',y;h)dx' < / e - h "™ {y - x )/e dx' < / e- x dx' < - , 

J x ^min JO ^min 

and the result follows. □ 



Proof of Lemma \4-8\ By the definition of p, the argument of the exponential defining </> is every- 
where nonpositive so <f> < 1 and I(p; h) < 2A1 . Moreover at the point p, the primitive of h is 
locally Lipschitz continuous, so there exists a sufficiently small r\ and a constant c > such that 

< — / h(x') dx' < -\x — pi, Ix — p| < 77. 
J„ c 



Then 

I(p;h)> 



[ e dx = 2ce f e x ' dx' > ce. 

J-n Jo 



-v 

□ 



Note that in Lemma l4.7l if /i m i n = 0, then the considered integral remains bounded as e vanishes. 



The following result Lemma 14.91 shows that, on any compact subset of the complement set 
[Amin, A max ] c , the mass of the linearized wave measures tends to zero. In the limit, all the mass of 
the wave measure <p(-;h) is concentrated on the interval [A m i n , A max ]. 

Lemma 4.9 (Behavior of linearized wave measures). For all h in the class L the function ip(-; h) 
satisfies the estimates 

( e -(z-A mln ) 2 rf min /2 £; -M <x< A min , 

< ip(x; h) < 0(l/e) < 1, x £ [A min , A max ], 

[ e -(^A max) 2 d min /2 % Amax <x< M _ 

Proof. For x > A max we have, with p defined through Lemma 14. 5[ and d and A given in Defini- 
tion S3 

h(y)dy = / d(y) {y ~ \{yj) dy + / d{y) (y - \{y)) dy 

"'A max J Pi 

> J* d(y) (y-A max )rfy> (x ~ Am 2 ax)2rfmi " , 
while a similar argument for x < A m ; n gives 

Hy)dy>J* d(y)(y-X min )dy> {x ~ Xm f 2dmin . 

Finally, the desired conclusion follows from definition of ip(-; h) and from Lemma FOl □ 

4.4. Wave coefficients. We rely on the earlier work by Joseph and LeFloch [35], in Lemma f4. 1 01 
Our new contribution is about the case of resonant wave coefficients treated in Lemma |4. 121 

Lemma 4.10 (Estimates of the wave coefficients). The linear coefficients Jj^,i defined in (|4.21|) 
the estimate 



(4-26) \Jj^{y)\ < 



0(e)(ipt(y)+^(y)), i + j, 

for all i,j = 1, . . . , N and y £ [—M,M]. Moreover, the quadratic wave coefficients defined in 
(14~22"1) satisfy 

(4-27) l^.k-xG/)! <C(<p*(y) + V](y)+<fk(y))- 
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Lemma 4.11. By choosing 5 small enough, for all i ^ j there exists positive constants C and D 
independent of e such that 



(4.28) 



< Ce- D '\ 

Lemma 4.12 (Resonant quadratic coefficients). Given rp G L 1 and i ^ j, then one has 

Jt^{y)\<omnM{y) + ti(y)l ^>o. 

Proof of Lemma \4-10\ The case i = j is obvious, so we only need to consider the case i ^ j ■ For 
definiteness we suppose that j > i, the proof for j < i being similar. 
First, using Lemmas 14.71 and 14.81 we get in the region y > c^: 



<pUv) 



dx 



v>j{y) / 4>{x,y\iH - m) dx 



< 



A, -A 



On the other hand, in the region y < Ci we have 



<pt(y) 



dx 



1 J Jy 



But, by Lemma 14.71 we have 



Vi(y) Y 4>{pj, Pi', Pj) <t>{pi,Ci; Pj - Pi 



(ci, x, pj - pi)dx < 



,x,p 3 - pi)dx. 



4i - A * 



and, by an easy computation, 

4>{p h c i ;p j -p H ) <e^-A i )|c i -«IA ) 
Ji/I, = 0(l/c). 
Using these observations we also get for all y < Ci 



V ^dx 



if*(x) 



with 



ft; 



-(Aj - AJ |cj - 



CAj — Pi) 2 



2(1 + 77) 

> — (Aj - A,) (Aj - A ( ) + 7 i 7 T ^y(A J - A,) 2 

When (5q tends to 0, (Aj — A.J remains bounded, (Aj — AJ vanishes, while (Aj — Aj) tends to 
Aj(0, 0) — Ai(0, 0) 0, so assuming <5o small enough, we can suppose each quantity /3jj is positive. 
The first desired result (|4.26[) therefore follows. 
Using the inequality \ab\ < (a 2 + b 2 )/2 we have 

\Fj, k -+i\ < ^(I^-mI + I^m-mD- 
So, we only need to consider the coefficients of the form Fjj^i, that is, 



4>{pi,y,Pi) 



(/)(p 3 ,x;p J )' 2 
(j>{pi,x;m) 



dx 
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To avoid any distinction between the cases y > Cj and y < Ci, the integral would be noted J, , 
(this is allowed by the positivity of the integrand) . Clearly, when j = i we have 

since 

4>(pi,x;fii) dx < Ii. 



So, we now suppose i ^ j, then 



_ <t>(pj,y;H) 
'I 



Gi,j(y) = Z^ll^ni I (j>(p i ,x;-iM)<l>(pj,x;iJ,j) 2 dx 
J ic>,y] 

<P(x,y,Pi - p,j)(f>(pj,%;iJ,j)dx 

,y] 

(/>(x,y;ni - p 3 )dx, 

,y] 



thanks to the judicious choice of pj that gives 4>(pj,x; pj) < 1. Note the important use of this L°° 
estimate that will default in the future. Finally we have 

Gi d {y)<T Jj^(y) 

which, together with (|4.26[) and Lemma 14.81 completes this proof . □ 
Proof of Lemma \4-H\ Fix y 6 [Aj, Aj] and recall that 



<pUv) <-exp(- / {\i(t)-t)dt 



If V > Pi then we find 



(Ai(t) -i)dt< / (A, - t) dt < i ((A, - Pl f - (A, - yf) 



< 



1 



((M-Arf-iM-yf) 



while, if y < pi, 
rv 

(Xi(t)-t)dt= I (t-Xi(t))dt< I (t-AJdt 



< \ ((Aj - Pi ) 2 (A, - yf) < ± ((A, - Aj) 2 - - yf) 



If i < j then for all e > 0, y < p% and 

rv j 



(Ai(t) -*)#<- ((A, - A,) 2 - (A, - A,.) 2 ) 



while if j < i then for all e > 0, y > pi and 



'(A. i (t)-t)^<i((A,-A l ) 2 -(A l -A,) 2 ) 



In all cases, we thus can write 



^^expQ-^-A 2 ,) 



where ti = Aj — Aj represents the width of the i-th wave and Ajj the gap between two different 
waves Ajj = min(|Aj — A_.|, \Aj — AJ). 
Moreover, we have 

< 2Mexp (-- f V (Xj(t) -t)dt) < 2Me (K '-^' i < 2Afe £ ? /e , 
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so finally, 



< 2MCe~ 1 e : 

<Pi(V) 



Choosing Sq small enough, we can assure the positivity of all quantities A? - — £j and the result 
follows. The multiplicative coefficient e _1 can be absorbed by reducing the exponential factor. □ 

Proof of Lemma \4-.12\ Using integration by part and denoting "J/ an anti-derivative of ip, (with so 
Halloo < HV'lli < we obtain 



dx \if*(x) 



^{y)^{y)-^{ymoi)^ ] 



<Pi(<H) 



dx 



dx \ip*(x) 



dx. 



The explicit formula for ip^ gives 



dx \ip*(x) 



d_ 

dx 



I /l C x | r x 

Y^v{- {\j(t)-t)dt~- (Xj(t)-t) 

1 <Pt(%) 
(A,(x)-A,(,))^ 



dt 



and, consequently, 

JUM <II*IIcoK(»)+^) 



|Ai-A 



3 Woo 



<pt(y) 



dx 



Thus, by Lemma 14.111 and Lemma 14.101 on binary terms for different wave families, 

\J 3 My)\<0(e) (^(y) + 

we get 

JUM\ < Mi (viM + <Pi(v)Ce- D ' e + 0(l)(^(y) + ^{y))) ■ 



□ 



Remark 4.13. The method just employed could be used to establish Lemma \4-.10\ directly, noting 
that H^JUl 1 = 1 an d Jj-Xi — Fj.k^i, and on the other hand, ||^||x,i = 1 o,nd J^i^ — Fj^^i- We 



deduce 



and (|4T77|) follows. 



\F j , k ^ i (y)\<o(i)(9*(y) + v*i(y)), 

\Fj,k-*i(y)\<0(l)(<p* k (y) + <pUv)), 
5. Construction of the entropy solution 



5.1. Correction vector for a given strength. Let Co(R) be the space of all continuous func- 
tions that decay to zero as |£| —> +oo, and define the following weighted sup-norm for 6 € [Cc^R)]^ 
as 

N 



Thus we search for 9 in the Banach space 

(5.1) E={e = (9 1 ,...,d N )e[C (M.)} N ■ 

For S > 0, consider the ball with radius S 

B s := {reR N : |r| < 6} 



< oo 
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and for r 6 B$ 

(5.2) |0 fc (OI < M V \t\ + |r| 2 + v\r\) ]T ^(£), fc = 1, . . . , Jv}, 

where A is a positive constant to be chosen later and v := Halloo. 

The set 5F is a closed bounded subset of E in the weighted norm ||-||. The quadratic quantity |r| 2 
is already present in Tzavaras |34] and comes from quadratic coefficients between the Ti^p* waves. 
Note however that the presence of the coupling wave ip (of unit total mass) , with its strength v, 
enforces the subset to contain correction waves that comes from coefficients associated with tp 
and the Tj(/J*, and of strength at most v\t\ relative to the <p* h . 

Now, we will define, for a given strength r the correction 8(t; •). Let define the map T that 
takes u £ f2 where 

n := {u £ C°([-M, M]), |K) - «£!!«, < <r} , r e B« 
and 8 £ 3 to the vector- valued function T(u, r, 0) whose components are given by (k = 1, . . . , N) 
T k (u,r,6)(0 

(5 ' 3) +V*(0 / ^^I^«ijfe(^)(^-(^)+^(x))(T i ^(a : )+^(x))dx 

+ <£*■(£) / —^r-ry2^k(x)(T l ip*(x)+8 l (x))^(x)dx. 

Lemma 5.1 (Contraction property). There exists positive constants A, r\, So and v such that for 
S < S : 

(1) T : £l x B$ x 3 3 is well defined. 

(2) There exists < a < 1 such that 

\\T(u,T,d) -T(u,t,6)\\ < a\\6-§\\, 8,8 £ J, 

and for any u £ fl, t £ Bg. Therefore T(u, t, ■) : $ J is a uniform contraction. 

(3) There exists a positive constant C , depending on fi but independent of 8, such that 

\\T(u,T,8)-T(u,r,d)\\ < C{ii + v + S)\t-t\, t,t£ B s 

and for any u £ Q, 8 £ J. 

We deduce from this lemma the following existence result of a correction 8(t: •). 

Proposition 5.2. Given u £ f2,r £ Bg. there exists a unique 0(t; •) £ J, i.e. in the class of 
functions satisfying 

(5.4) \6 k (r;-)\ < A( V \t\ + \t\ 2 + v\r\)J2<Ph, |r| < S, k = 1, . . . , N, 

h 

solution of the fixed point equation T(u, T, 6) — 8. Moreover, there exists a constant C independent 
of 5 such that 

(5.5) \d k (T;-)-8 k (T;-)\<C( V + v + 5)\T-T\J2vh r,f £ B s . 

h 

Proof of Lemma \5.1\ The main difficulty is handling the coupling wave wip. First of all, we show 
T keeps the subset J stable, using the definitions (|5.2I) . (I5.3[) . and the different definitions of the 
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coefficients (H^Tj) . (|3^2l and we get 

|T fc (u,r,0)(O| 

< vMooYl "WOfVl + A ( V \ T \ + \r\ 2 + u\r\) ) 
+ || K |U £ F^ fc (0(|r| 2 +2|r|A(r?|r| + |t| 2 + v\t\) + A 2 (rj\r\ + |r| 2 + V \r\f) 
+ Hoc J2 J tk(0(\r\ + A ( V \r\ + |r| 2 + u\t\) ) . 

i 

Thus, by using Lemmas 14.101 and 14.121 we get 

\T k (u,r,6)(0\ 

< vlWWoo C,N 2 (|r| + A ( V \r\ + |r| 2 + u\r\)) £ 

+ ||«l|oo C 2 JV 4 (|t| 2 + 2A|r| (771-7-1 + |r| 2 + i/|r|) + A 2 (^|t| + |r| 2 + u\r\) 2 ) ]T 

ft 

+ Hoc C7 3 iV 2 (|r| + ,4 (771-7-1 + |r| 2 + u\r\)) ]T ^(0 

h 

<C (I + A(ri + v + 5)f (rj\r\ + |r| 2 + V \ T \) £ c^(0, 

/i 

where C is a constant depending only on AT the dimension of the space, on ||7r||oo, ||«|| oo an d of 
the constants Ci, C 2 , C3. A necessary condition to get the stability of the subset J by T(u, r, •) is 
also 

C(l + A(r) + v + S)) 2 < A 
A way to get this inequality is for example, fixing A — AC, to choose 77, v and 5 together such 
that rj + v + S < 1/4C. 

Now, T is an uniform contraction relative to the variable (?e J, since (from similar arguments) 
\T k (u,T,6)(0-T k (u,T,§)(0\ 

<^Hoc]Tll^llE J ^) 

3 i 

+ ||«||oo ^(2|r| \\0 6\\Y,Fu(0 + 2A ( V \r\ + |r| 2 + v\r\) \\d 0|| 5> ItB -*(£)) 

ij I Im 

+ h\\ooJ2w e - § \\J2 J ^k(0- 

3 i 

In view of Lemmas 14.101 and 14.121 we obtain 

\T k (u,r,em-T k (u,r,0)(O\ 

h 

+ Moo c 2 n± (\t\ + a (t)\t\ + \t\ 2 + v\t\)) \\e -ewY, <pUZ) 

h 

+ IHU C 3 N 2 \\6 - 6\\Y,vl{£) 

h 

< C(ri + v + S)(l + AS) \\9 - §\\ 

h 

Finally, we obtain the item (|2|) of Proposition 15.11 with a — C(rj + v + S)(l + AS) by choosing 
rj + v + 5 sufficiently small to assure that a < 1 . 



22 



BENJAMIN BOUTIN, FREDERIC COQUEL, AND PHILIPPE G. LEFLOCH 



Finally, we have to check the Lipschitz continuity of T in the variable r. 
\T k (u,T,8)(0-T k (u,T,6)(0\ 

<vh\\oo\T-f\J2 Ji ^ k ® 



+ N|oo53 (infj -TjTilFij-n,® +2|r -f| A5(r] + S + u)^2Fu^kiO 

ij I 

+ Woo|r-f|X;^*(0- 

i 

By using Lemmas 14.101 and 14.121 we get 

\T k (u,T,6)(O-T k (u,t,0)(O\ 

< vMoodNlr - f | ^ ^(0 + \\k\\ 0O Q 1 N 9 (2S\t - f | 

h 

+ 2ANS\t - f |) ^ <p* h {£) + \WWnCaN\T - f 1 53 vUO 

h h 

<C(r) + v + 5)\T-t\ 53^(0, 

h 

and the item ^ of Proposition 15.11 follows. □ 



Proof of Proposition \5.Sl The inequality (|5.4j) is a direct consequence of the contraction mapping 
theorem, that previous proposition ensures to apply (a < 1). Let u 6 f2 be given, and r, f G -B,5, 
then 

0(r) - 0(f) = T(«, r, 0(r)) - T(u, f , 0(f)) 

= (T(u, r, 0(t)) - T(u, r, 0(f))) + (T(u, r, 0(f)) - T(u, f , 0(f))) , 

thus 

\\0(t) - 9(f)\\ < \\T(u, t, 0(r)) - T(u, r, 0(f ))|| + ||T(u, t, 0(f)) - T(u, f , 0(f ))|| 
< a||0(r) -0(f)|| + C(t) + v + 5)\t -t\. 

Hence, 

||0(T)-0(f)||<-^-(77 + ^ + <5)|T-f|, 

1 — a 

and (15.51) ensues. □ 



5.2. Strength vector for given Riemann data. Fix a left-state vector uz G and u G Q = 
{u G C°([—M,M]), \\u(-) — wl||oo < c}- Being given r G _B^, we previously constructed a unique 
0(t, ■) £ J such that T(u, r, 0) = 0. The question is now to link the vector t e Bj to the boundary 
data ul, ur. Consider the following operator: 

rM 

S(t) := Ul + A ( W , U )- 1 53 / [r W UO + O k (r,0}r k (u(0,v(0,m- 

k J - M 

Lemma 5.3. There exist constants 5, r > such that the operator P : B r (uL) x Q x Bg — > P>s 
defined by 

l-M 

P(u r ,u,t) = Aq(u,v)(u r -u l ) -53 / ®k(T,£)r k (u(£),v(£),£)d£. 

k J - M 

satisfies the contraction property (for some < a < 1 ) 

\P(ur, u, t) — P(ur, u,t)\ < a\r — t |, r, f G Bg 
for any ur G B r {ui),u G Cl. 

Proposition 5.4. Given u E Q, there exist positive constants r and S such that the following 
holds: 

(1) For all ur G B t {ul) there exists a unique solution of the equation S(t) = Ur with t G Bg. 
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(2) For each u £ fl and e > the inverse map S 1 : B r (uif) — > Bg is well defined and satisfies 
(5.6) \S~ 1 (ur)\ < i\ur - u L \, 

where 7 is a constant depending on q but independent of u £ Q and e. 
Proof of Lemma \5.3\ Letting ur £ B r (uif), u £ Cl and r £ Bs, one has 

\P(u R ,u,-)\ < \\A \\\u r -u l \ + N 2 RA(ti\t\ + \t\ 2 + u\t\) 
< (\\A \\r + RN 2 A(ri5 + S 2 + vS)) . 
Hence, one has the inclusion P(ur,u, Bs) C Bs provided 

(\\A \\r + RN 2 A{ V 8 + S 2 + v5)) < 8, 
that is by choosing r, T], 8 and v such that 

RN 2 A(?] + 8 + v)8 < 8/2, 
\\M\0r < S/2, 

that is to say 

ri + 8 + v < 1/2RN 2 A, 
r< 8/2\\A \\. 

Given r and f in Bs we have 

,-M 



rM 

\P(u R ,u,r)-P(u R ,u,t)\<J2 / \9k{Z,T)-6 k {Z,T)\ \r k (u(0,v(0,0\d£ 



-M 

pM 

< RNC( V + 8 + v)\t -t\Y, <pU® # 

J-M 



k 

t2, 



< N z C{r 1 + 8 + v)\T-T\. 
Provided 

(5.7) a := N 2 C(r] + 8 + v) < 1, 

the map P(ur,u, •) is a uniform contraction on Bs. □ 

Proof of Proposition \5~^\ Let ml be fixed. The equation S(t) = ur takes the form 
A (u,v)(u R - u L ) 

/M r-M 
-M k J-M 

in other words r solves the equation 

pM 

(5.8) A (u,v)(ur-u l )=C(u,v) t + Y, / h(T,0?k(u(0,v(0,Odt, 

fc J - M 

where C(u, v) is the matrix whose fc-th columns is given by 

<pm?k(u(OMO,zn, k=i,...,N. 

M 

This matrix has the important property it is invertible for any u £ f2 and the inverse matrix 
C(u,w) _1 is uniformly bounded (cf |34j ) 

(5.9) IC'iu^)-^ < p, u£Cl. 

In order to solve the equation S(t) — ur, observe that solutions of (|5.8|l are also fixed points of 
the map r i-> C(u, v) P(ur, u, t), whose existence are ensured by Lemma I5T31 As a consequence, 
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given Ur G B t (ul), there exists a unique fixed point r of P(ur,u, •) in the ball B$. Moreover it 
also satisfies 

\t\ < \A Q (u,v)C(u,v)~ 1 \ \u R - u L \ 



/Ivl 
\0k^,r)\\r k (u(O,v(O^M, 
-M 



-M 

thus 

\r\ < \\A \\P\ur -u L \+ /3RAN 2 ( V \t\ + |r| 2 + v\ T \) 
< \\Ao\\P\u r -u l \ + 1/2\t\. 
Thus, |r| < 2\\A \\/3\u R - u L \, which finally implies ([BTB]) . □ 

5.3. Riemann problem. We search for a solution of (|4.14l) under the form (|4.15l) satisfying the 
boundary conditions (|4.12[) and where v is known by (|4. 13|) . 

Theorem 5.5 (Uniform estimates and existence result). There exists a solution u f (^Qof the 
problem (|4. - (|4. 14p satisfying, for some constant K independent of e, 

TV(u*)<K\u R -u L \, 

(5 - 10) e\u\\<K. 

After extracting a subsequence if necessary, this result provides us with a solution with bounded 
total variation. 

Proof, q > is choosen so that conditions of eigenvalue separation are fulfilled on f2. Fix ul and 
u G Q. For e fixed, we construct z as 

z(0 = u L + A (u,v)- 1 / X>i$(O + 0,(r;O)fi-MOMC),CK 

= 1 

by following steps: 

(1) Each if* is constructed as the fundamental wave measure from (|4.18j) . recalling 

(2) For each r small enough we can get, through Proposition 15.21 a correction 0(t, •) so that 
dj — Tj ip* + 6j is solution of (|4.16a[) 

a' — — — dj = r]Lj(u,v, •) + (u, v, •) + ^(u, v, ■). 

J e 

(3) The vector of strength r is then chosen, through Proposition 15. 4[ as a solution of S(t) = 
uji G B r (uL). This way the solution u of 

A (u,v)u^) = a 3 -(C) ?,•(«,«,£) 
j 

satisfying u(— Af) = satisfies moreover u{M) = Ur. 

These steps allow us to construct an operator 7 : f2 — > E, u i-> z, and T(w) = 2£ft. We only need 
to get a fixed point result on 7 to get the solution u of the whole problem, and then sufficiently 
strong estimates to ensures existence of the limit as e tends to 0. □ 

Lemma 5.6. The function v e converges toward the sign function (denoted by sgn) and, more 
precisely, for all c > 

(5.11) \\v £ - sgn\\ Loa{m _ c c]) = o(e). 

Proof. Indeed the formula (|4.13|) implies v e takes the form 

J_ M e ^dx 



COUPLING TECHNIQUES FOR NONLINEAR HYPERBOLIC EQUATIONS 



25 



Fix £ > c > 0, so that 



1| < 2-^ 



rM _ , 

J_Af e 



< 2- 



£ p/2 



For £ < — c < 0, by the same procedure, we get 



C 

< — 7^e~^ . 



p/2 



\v e (x) + l\ < ^e"^, 



the lemma is therefore proved. 



□ 



Theorem 5.7 (Convergence to an entropy solution). The sequence u e converges pointwise toward 
u G BV , satisfying 



(5.12) 



-e^7-(«) + 4/-(7-(«))=Q, mt<0, 
< ^7+(«) + ^/+(7+(«)) = 0, m£>0. 



Let ?7± = ?7± (it) G &e two entropy functions compatible with the viscosity matrix in the sense 
that 

V 2 r] ± (u)B ±(u) > 0, it e ft. 
Then following entropy inequalities are satisfied, 



(5.13) 



-C ^»7+(7+(«)) + ^?+(7+(«)) < 0, m(>0. 



Proof. Let $ be a neighborhood of 0, the solution u c of (|4. is obtained as a smooth function 
so that, under the consistancy hypothesis (|2.10[) , we have outside $ 

Ao(u>0) ^ = ! 7±(u£) ' 
^iK,-°)^ = |/ ± (7 ± (- e )). 

Let (/) £ C^°(K_ \ $) be a test-function with a compact support included in K_ \ then (|4.1ip 
implies 



-£ (AoK,«°)-A o ( u £ ,« £ )) u |0^ 

+ / (A 1 {u € ) v°)-A 1 {u e ,v e ))u\tt>d^ + f e(B (u e ,v e )u\) <j> d£. 
Moreover, thanks to Lemma 15.61 and using Lipschitz continuity properties, we have 

-£ (Aq{< v°) - A (u e ,v e )) u\cj> d£ < O (e)Lip(A )||^||ooTV r ( M e ), 
(A x {u e , v°) - Ai(u £ , v e )) u\cj> d£ < (6)Lip(A 1 )||</)|| oo TF( W e ), 



and 



(B (u 6 ,u>£)^d£ 



< e||Bo||oo||^||ooT^(« e )- 
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Thus, as e tends to 0, we get the weak formulation for the limit u 

H ^7-(«) + ^/-(7-(«)))vde = 0. 



By a similar method we get for ip G C° 



H ^7+(tt) + ^/+(7+(«)))vde = 0. 
Entropy inequalities are obtained by first using the consistency hypothesis (|2.10p . that give outside 

V7?±( U £ ) • A {u\v°)u\ = ^r?±(7±(n £ )), 

V77±K)-AiK,«°) W | = ^ ± (7 ± ( U e )). 

Let (/> G C^°(K_ \ $) be a non-negative test function with a compact support included in R_ \ <I>, 
then (|4.1ip implies 

H ^-(7-(« e )) + ^?-(7-(« e )))^de 



+ / V»7_(u £ )-(Ai(u e ,« o )-Ai(^» e ))u!0d£ 

+ / eV V -(u e )-(B (u e ,v e )ul) <j>d£. 
Jm. 5 

With similar arguments as previously, the first and the second terms of right hand side tends to 
as e tends to 0. Moreover, after reporting the ^-derivative on (/>V?7_(i/), the last term 

^V)j_(u e ).flo(« e s »>|de- / #VV(« £ ) • B {u\v a )\u\\ 2 d£ 



+ / e<jN 2 r,-{u*) ■ {B {u\v Q ) - B {u*,v e ))\u\\ 2 d$ 



satisfies the estimates 



and 



/ e0 5 V77_(u e ) • B (u e ,v e )ul d£ < KeTV{u e ) 



e0VV(" e ) ' (B (u e ,v°) -B (u e ,v'))\ul\ 2 d^ 

< jaipCSoJIIvVllooll^Urv^iK-^iu-^*). 

However the quantity J R e</)V 2 ry_ (u c ) • Bo(u e , u°)|w|| 2 d£ is not guaranteed to vanish as e tends to 
0, but it converges toward a positive value under the hypothesis V 2 rj±Bo±(u) > 0. The following 
weak formulation of the entropy inequality on R_ follows: 

K ^T7_( 7 -(«)) + ^g_(7-(«)))0 ^ < 0. 
Similar arguments lead to the entropy inequality on R + . □ 
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